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Abstract 

Numerical results for Anderson transition are critically discussed. A simple 
procedure to deal with corrections to scaling is suggested. With real un- 
certainties taken into account, the raw data are in agreement with a value 
v — 1 for the critical index of the correlation length in three dimensions. 



Critical indices s and v for the conductivity o and the localization radius £ of wave 
functions near the Anderson transition [Q are determined by relations 

a ~ r s , £ ~ r~ v , (1) 

where r is the dimensionless distance to a critical point. On the modern level, the first 
estimates of s and v were given by a scaling theory of localization 0, s = 1 and v = 1/e for 
the space dimensionality d = 2 + e. They suggest that s ~ v ~ 1 for d = 3. A little later, 
a self-consistent theory of Vollhardt and Wolfle || was developed: it gave for arbitrary d 



1 , d > 2 



2 < d < 4 



d-2 

(2) 

\ ' d>4 

These values of indices were in a good agreement with all existing information and suspicion 
arised that they are exact . Later the present author has suggested a phenomenological 
scheme, based on the symmetry considerations, that gives (2) without model approxima- 
tions. 

Results based on a nonlinear a-model M agree with (2) on a three-loop level, but the 
correction of order e 3 appears to be finite |J and shifts s and v to much lower values, 
strikingly spoiling agreement with experimental and numerical results. One should have 
in mind, however, that correspondence of a nonlinear cr-model with the initial disordered 
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system is approximate and valid only for d = 2 + e with small e. So it is not surprising if a 
difference arises in a certain order in e. On the other hand, the high-gradient catastrophe 
|]] makes the essential modification of the cr-model approach to be inevitable [[TTJ . 

A value s = 1 is in agreement with a large number of experiments [11] but this agreement 
can be doubted on the ground that the interaction effects are essential in real physical 
systems. However, attempts to include the interaction in the Vollhardt and Wolfle scheme 
[0 show that a result (2) can preserve in the systems with interaction. Convincing evidence 
for s = 1 was obtained recently in the experiments with the nondegenerate electron gas 
[|j~3|l where the interaction effects are negligible. 

Early numerical results were in a reasonable agreement with (2) [y = 1.2 ±0.3 |TJ|], v = 

but later a tendency to the larger values arised: v 
1.54 ±0.08 



0.9±0.3, v = 1.4±0.2 [p 
= 1.50 ± 0.15 pL 



IE 

v = 

indices. 



v 



1.45 ± 0.08 M 



v 



1.35±0.15 



1.4 ± 0.15 



1.58 ± 0.02 [f2T|. So large values contradicts to all other information^ on the critical 



It means, in our opinion, that something is wrong with numerical results. 

1. To understand it, let us consider <i-dependence of the numerical value for u, which 



is given by an empirical formula [23] 



0.8 



v 



d 



±0.5 



(3) 



This formula summarizes results for d = 3, d = 4 and several noninteger dimensionalities 
which are realized in the fractal structures. From the theoretical viewpoint, the formula 
(3) is entirely unsatisfactory: 

(a) It gives v = 0.8/e for d — 2 ± e, while it should be v = 1/e in any variant of 
one-parameter scaling: it is a consequence of the fact that the Gell-Mann - Low function 
has a behavior (3(g) = (d — 2) ± A/g for large g, the latter property being proved by a 
diagrammatic analysis. Of course, one can doubt in validity of one-parameter scaling, but 
in such case the whole procedure of the data treatment should be dismissed as entirely 
based on it. 

(b) Formula (3) takes into account that v = 1/2 for high dimensions [|4] but suggests 
the infinite value for the upper critical dimensionality d c2 . Such hypothesis arised in the 
formalism of a-models, which do not exhibit any special dimension except for d = 2 p5 



Once again, it is a consequence of the fact that the a-model approach is approximate and 
can be justified only for low dimensions. 



1 In a recent communication 22 A. Kawabata claims that a value v = 1.58 can be obtained, if self- 
consistency of q-dependence of the diffusion coefficient D with L-dependence of conductivity is required. In 
fact, a solution with v = 1 is already self-consistent in this sense: it has no g-dependence, but L-dependence 
of conductivity follows from the temporal dispersion of D ||. The general analysis shows Q that only a 
solution with negligible q-dependence is self-consistent. 

2 Efetov's derivation ]25| ] of the cr-model is formally valid for arbitrary dimensions but includes the artifi- 
cial construction of weakly-coupled metallic granules. Such construction and related with it approximations 
can have unpredictable influence on results. 
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In the exact field theory formulation, the theory of disordered systems is equivalent to 
the ip A theory with a "wrong" sign of interacton [26, 27]. The latter is renormalizable for 



d < 4 and nonrenormalizable for d > 4 [2B\ []. For d < 4 all of the physics is determined 
by small momenta or large distances, in accordance with the expected scale invariance. 
For d > 4 the atomic scale cannot be excluded from results and no scale invariance is 
possible. So the upper critical dimension is surely four and it can be seen from the different 
viewpoints |3(J. So we should expect v = 1/2 for d = 4, but it is in conflict with numerical 



value v = 1.1 ± 0.2 ||31| , which can be doubted on these grounds. By continuity, the 



analogous systematic error can be expected for d = 3. 

2. The most probable reason for a systematic error is the existence of corrections to 
scaling. All numerical results are based on one-parameter scaling and indeed there are no 
serious doubt in it for low dimensions. However, there are no indications of the existence 
of the upper critical dimensionality in the framework of one-parameter scaling. The only 
possibility to resolve this controversy is to suggest that one of irrelevant parameters becomes 
relevant at d = 4. The corresponding scenario was developed in [^] and resulted in the 



scale dependence of the conductance which is in agreement with the Vollhardt and Woffle 
scheme. There is two-parameter scaling for d > 4 and a large systematic error for d = 4 
is rather natural. For d — 4 — e, two-parameter scaling takes place for intermediate scales 
and the region of its applicability becomes smaller with increase of e. As a consequence, 
disagreement of numerical results with (2) is less for lower dimensions. 



On the model level, corrections to scaling were analysed in |p3[ . With some assumptions, 
it appears to be possible to agree the raw data with a value v — 1. 

As it was claimed in with corrections to scaling taken into account, the accuracy 
of results for v reaches the level ~ 1%. This statement is extremely doubtful. Authors of 
fJIjj carry out nonlinear fitting with 10-12 parameters and present only results but not the 



procedure used. But in such situation, \ 2 nas a great number of different minima (as in 
a spin glass), and the true one is not necessary the deepest: there are no real possibility 



to systematize them. Probably, authors of [21] have in mind the minimum, nearest to the 



fit without scaling corrections. But a problem is, that such fit may be not a good zero 
approximation. 

3. We can suggest a simple procedure to deal with corrections to scaling. As a starting 
point, we suppose existence of the abstract renormalization group, which is determined by 
the operator R s : it corresponds to a change of the length scale by a factor s and transfers 
one point /i of the parameter space into another point //, // = R s fi ||26|| . If fi* is a fixed 
point, /i* = R s fi*, then we have for small 5fi = fi — /i* 

5fi = fi - fi* = ^Ai{r)ei and R s 5/i = ^ A i (T)s Vi e i , (4) 

i i 

where and s Vi are eigenvectors and eigenvalues of the operator R s , considered as linear 

3 Problems of the "wrong" sign of interaction and the replica limit n — > are inessential in the 
framework of perturbation theory. In fact, renormalizability can be investigated directly in the framework 
of the "impurity" diagrammatic approach p9[, without transition to the effective field theory. 
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in the vicinity of \x* ||26|| . If \x is represented by a set of the coefficients Ai(r), then we can 
write for any quantity Q with a zero scaling dimensionality: 

Q(ji) = F{A 1 (t)s*, A 2 (r) S y\ A 3 (r)^, ...} , (5) 

where y\ > y 2 > y% > . . . and y\ is related with the critical index v of the correlation 
length, yi = 1/v |36]. If we expand Aj(r) near the critical point, Aj(r) = -B^r, we have 
after redefinition of F 

Q(r,L) = F{r(L/a)y\ r(L/a)^, r(L/a) w , . . .} , (6) 

where a length scale L is s times greater than an atomic scale a. For r(L/a) yi <C 1 all 
arguments of F are small and 

Q(r, L) w F {0, 0, 0, . . .} + driL/a)* 1 + C 2 r(L/a) y2 + C 3 r(L/a) y3 + ... 

= F {0, 0, 0, . . .} + rf {L/ a) , r{LjaT « 1 • (7) 

In the case of one relevant parameter we have y\ > 0, y 2 , y3, ■ ■ ■ < and for large L/a 
only the first argument in (6) produces essential changes in F, while influence of other 
parameters are uniformly bounded by a quantity of order r: 

Q(r, L) w F {r(L/a) w , 0, 0, . . .} = G {r(L/a) w } , r(L/a)^ > 1 . (8) 

Now we can unify (7) and (8) into one expression: 

Q(r,L) = G{rf(L/a)} . (9) 

A function / has a power-law asymptotics for large L 

f(L/a) ~ (L/a) 1/u , L/a -> oo (10) 

and an arbitrary behavior for L/a ~ 1, when all terms C{t(L/ ' d) Vi in (7) are of the same 
order. Nevertheless, dependence of In /(re) on lnx is slow for a wide class of smooth 
functions and in a restricted interval it can be linearised. So 

f(L/a) w const (L/a) 1 ^ (11) 

and we have a relation 

Q(r,L) = GfrCL/a) 1 /"-"} = G{i{r)/L} , e(r) ~ ar"^ , (12) 

which is suggested in all papers involved in the numerical business. At L — » oo index u e ff 
tends to u, but at finite L it may be essentially different. 

Expression (9) is more general than (12) (only r < 1 was suggested) and a simple 
procedure of the data treatment can be based on it. The most extensive information on 

4 We assume that a function F and its derivatives are bounded. 



4 



the function G can be obtained for a maximum value of L, which is available for us. If we 
accept f(L max /a) = 1, fixing the scale of the argument of G, we have 



G{r} = Q(r,L max ) . (13) 

If Q(t, L) is known for a set of the scales Lj, the formula (9) determines f(x) in the points 
Li/ a: we should fit the scale Tj to a relation G{r/rj} « Q{r,Li) and put /(Lj/a) = l/Vj 
after that. Such procedure makes it possible to extract the true function f(L/a), not 
assuming a power- law dependence for it. Its behavior for L ~ L ma:r gives the most reliable 
estimate of u, that can be obtained from the existing data. 

4. As an illustration, let us discuss the interpretation of the largest scale data, existing 
for the Anderson transition. The largest system size L = 100 was used in |2{J (we accept 
a = 1 in what follows), but the suitable data are available only for L < 28. The most 
detailed of them are presented in Fig. 4 of PD[ , where a scaling quantity A is plotted 
against the amplitude of disorder W for L = 28, 12, 6; we consider A as Q(t, L) with 
r = (W — W c )/W c , where W c is a critical disorder. 

Theoretically, the simplest procedure (in fact, it was used in p0[ ) can be based on the 
formula (7), which determined f(L/a) as a slope of a linear dependence Q on r at r — > 0. 
In practice, for sufficiantly small r one should take a difference of two close quantities 
with a great loss of accuracy, while for a moderately small r it is difficult to control the 
validity of (7). As one can see in Fig. 4 of |[2Q|| , dependence Q on r is only approximately 
close to linear, but in fact it is an essentially broken line. To demonstrate a situation, we 
give in a Table an average slope of dependence Q on r and its fluctuations in the interval 
16 < W < 17, which corresponds to a condition r(L/a) 1 ^ u < 1. 



Table 

Slope of dependence Q on r (arbitrary units) 





L = 28 


L = 12 


L = 6 


Average value 


0.30 


0.16 


0.10 


Least value for 16 < W < 17 


0.20 


0.10 


0.04 


Largest value for 16 < W < 17 


0.42 


0.25 


0.12 



With the use of the average slope, we indeed obtain v ~ 1.4, as it was reported in 
2~0|j . With real uncertainties taken into account, we can have any value of v in the interval 



0.7 3.0. Authors of |2(J give essentially smaller error, relying on the averaging procedure. 
But one can see by eye in Fig 4 of ||20|| , that any value, given in Table 1, has a reasonable 
probability to be a true slope at r = 0. Statistical treatment cannot improve this situation, 
but can give only an illusion of doing it. 

More stable results can be obtained with the use of a formula (9). If we accept /(28) = 1 
and take Q(t, 28) as G(t), we can find a scaling factor r, for each experimental point for 
L — 12 and L = 6. For a function f(x) we have 



/(12) = 0.33 4-0.65, /(6) = 0.14 4-0.40 



(14) 
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and after averaging 

(/(12)) = 0.52, (/(6)) = 0.28. (15) 

One can see in Fig. 1, that the most probable value is v — 1.25, and it is essentially shifted^ 
in comparison with v = 1.4 given in J20J. With the total uncertainty taken into account, 
we have v = 0.8 -j- 1.7 and a value v — 1 is admissible. 

Uncertainties of results can be formally diminished, if we interpret an error in the 
mean square deviation or another sense. However, it requires some statistical hypothesis 
concerning a distribution of errors, which hardly can be justified in view of their partially 
systematic nature. 

5. An effective index v e ^ in (11,12) tends to v in the limit L — > oo. For the moderate 
L, a good accuracy of approximations (11,12) does not mean, in general, that u e ff is really 
close to v. A well known example is the second order transition with a small Ginzburg 
number, where a new scale £ 3> a arises: the true critical indices are observed only for 
L <: £o> while for L < £ the relations (11,12) are valid with a mean-field value of v e ff- In 
fact, convergence u e ff to v can be slow even in the absence of such new scale. Evidently, 
there are no rigorous procedure to control a difference v e ff — v. 

All numerical estimates of v are systematically greater than unity and, beyond any 
doubt, it reflects a realistic behavior of f(L/a) in the corresponding range of scales. A 
function f(L/a) is not universal and depends on the choice of a model. One should have 
in mind, however, that one and the same model is investigated in all papers: it is the 
Anderson model with a transition in the center of the band. Differences in the form of 
distribution of cite energies (f.e. box or Gaussian) are of minor importance: a critical value 
of W is large, and near the band center this distribution looks as practically uniform. To 
have a real estimate of a systematic error, one should shift to a band edge and use the 
different procedures to reach a critical point (f.e. by change of energy or disorder). The 



only attepmpt of this kind was undertaken in [15] and resulted in essential decrease of v 
with increase of uncertainties. It is desirable to repeat such attempt on the higher level of 
accuracy. 

In conclusion, we have demonstrated, that even the most detailed and largest scale data 
can be reinterpreted with an essential shift of results. As for the systems of smaller size, the 
results can be doubted simply on the ground, that the function f(L/a) does not reach its 
asymptotic regime (10). As for a systematic error, nobody can say nothing definite about 
itfl 

5 The origin of this systematic shift is rather simple. The function Q(t) is linear for small r but grows 
more slowly, something like |r| a signr with a < 1, for larger r. If the latter dependence is approximated by 
linear one, using values for t = tq and r = —To, the slope ~ t§ is obtained. After a scale transformation 
r — > t{L/o) 1 / v this slope has dependence [L/a) a / v ', which is interpreted as (L/a) l ^ Ue ff with v e ft =vja> 
v. Such effect is present in the most of papers. 

6 The most extremistic values of v [|l8|, are obtained in the case of small L and high accuracy. It 
is exactly the situation, when a statisical error is small, while a systematic one is large. Surely, the given 
accuracy reflects only the first of them. 
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